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Abstract 

/ntroduction:T\r\e natural history of HIV-1 infection and its progression towards AIDS vary considerably among individuals. 
Host genetic factors may be one of the possible reasons for variable HIV-1 disease progression. Single nucleotide 
polymorphisms (SNPs) in the promoter region of TNF-sx gene can influence its production. The aim of the present study was 
to determine the association of functional TNF-a SNPs and its associated parameters related to apoptosis that may influence 
the rate of HIV-1 disease progression. 

Methods: Therapy naive, 100 HIV slow progressors (SPs), 100 HIV fast progressors (FPs), 50 HIV exposed but seronegative 
individuals (ESNs) and 260 healthy controls from same ethnic origin were recruited. Genotyping of TNF-oi variants {— 863C/A, 
-308G/A and -238G/A) was done using PCR-RFLP. CD4 counts were determined by flow cytometry. Plasma viral load was 
estimated by COBAS AlVlPLICOR HIV-1 monitor test. Plasma TNF-a concentration was estimated by Human CBA Thl/Th2 
cytokine kit. The lymphocyte mitochondrial membrane potential was measured by JC-1 dye by flow cytometry. 

Results: Genotype and allele frequency of TNF-a -238G/A and -863C/A was not significantly different in HIV-1-infected 
patients when compared to controls, while that of TNF-a -308G/A variant (high TNF-a producer) was significantly higher in 
FPs compared to SPs (p<0.01, OR = 3.43). Haplotype analyses also showed that carriers of high TNF-a producing haplotype 
CAG was significantly more common among FPs compared to SPs (p<0.01, OR = 3). The circulating TNF-a levels in blood 
also correlated well with genotypes. The lymphocyte mitochondrial membrane potential of FPs having CAG haplotype was 
significantly low as compared to wild type (CGG) haplotype {417±22 vs 571 ±28, p<0.01). 

Conclusion: High producer haplotype, CAG of TNF-a gene associates with enhanced apoptosis of lymphocytes in HIV-1 
infected individuals, hence faster progression to AIDS. However, further functional studies are needed to confirm this 
association and this knowledge may help clinicians to better understand the disease outcome. 
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Introduction 

The natural history of HIV-1 infection and its progression 
towards AIDS vary considerably among individuals. Polymor- 
phism of the 'immune response genes' including chemokine 
receptor/ligand system may be one of the possible reasons for this 
variable rate of disease progression [1]. Single nucleotide 
polymorphisms (SNPs) in the promoter region of cytokine genes 
may affect the quality of the immune response against a pathogen 
[2]. The impact of such SNPs on the course of HIV-1 infection has 
not been seriously investigated in Indian population. The tumor 



necrosis factor-alpha (TNF-a) is a pleiotropic cytokine that acts as 
an immune and inflammatory mediator, its gene located in the 
class III region of the major histocompatibility complex (MHC) on 
chromosome 6. Although its production is tightly regulated in vivo, 
the genotypic differences may influence transcriptional regulation 
[3]. The excessive TNF-a production is known to be associated 
with various autoimmune, inflammatory and infectious patholo- 
gies [4] [5]. Investigations suggest that TNF-a is involved in the 
pathogenesis of HIV-l infection, since it is overproduced by 
infected individuals [6] [7]. High concentration of TNF-a may 
influence HIV replication via clonal expansion of infected T 
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lymphocytes. In addition, TNF-a induces rapid apoptosis of 
bystander T helper cells, thereby accelerating the rate of T-ceU 
depletion in HIV disease [8]. Several functional SNPs have been 
identified within the T}fF-a gene cluster, the most widely studied 
of which are: G>A transition at position -238, G>A transition at 
position -308, and OA transversion at position -863 [9,10,11]. 
Few case-control studies have assessed the influence of TNF-a. 
gene-polymorphism on HIV-l disease progression in infected 
patients [12] [11]. While some of these studies have reported an 
association with disease progression [13] [14], others have failed to 
find such an association [15]. There is a lack of comparative data 
on the influence of TNF-a. SNPs on various cUnical cohorts of 
HIV-l infections like fast progressors, elite controllers, long-term 
non-progressors and exposed but seronegative individuals. 

The aim of the present study was to determine the associations, 
if any, of potentially functional TNF-cc gene SNPs, both 
individually and at the haplotype level, with the rate of HIV 
disease progression or risk of HIV infection in individuals from 
North India. In addition, we have attempted to explore whether 
any of these polymorphisms are associated with parameters related 
to apoptosis that may influence the rate of decline of CD4 cells. 

Materials and Methods 

1.1 Ethical Statement 

The study was approved by the Institutional Ethics Committee 
(lEC) of PGIMER, Chandigarh, India and an informed consent 
was obtained from all the patients and healthy controls before 
obtaining the blood specimen. 

1.2 Study Design 

Patients in different study groups were recruited from ICTC 
(Intergraded Counseling and Testing Center) in the Department of 
Immunopathology, PGIMER, Chandigarh, India. The study was 
carried out on HIV seropositive cohorts from North India that 
included: Group-I (Slow Progressors; SP, n= 100): therapy naive 
individuals who maintained CD4 count of >400 cells/ |xL for past 
3 years or more (slope of decline<-32 cells/|J.L/Year); Group II 
(Fast Progressors; FP, n— 100): individuals with fast depletion of 
CD4 count to <200 ceUs/|a,L within one year (slope of decline> 
-200 cells/(iL/Year). Group III (Exposed sero-negative; ESN, 
n — 50): individuals who were exposed to HIV but remained 
seronegative. Inclusion criteria in this group was: the HIV 
seronegative partners of discordant couples, having unprotected 
sexual intercourse for past one year with a minimum of 4 to 5 
exposures every month. The recruitment was based on a detailed 
history of each couple by the counselor to ascertain the 
unprotected exposure for a specified period of time and the HIV 
negative status based on serology as well as HIV-l DNA PGR. 
Group IV (Healthy controls; HCs, n = 260): HIV seronegative 
healthy volunteers from the same geographical region. HIV 
patients on anti-retroviral therapy (ART) and HIV patients with 
TB and other chronic co-infections (HCV, HBV) were excluded 
from the study. 

1 .3 DNA Extraction 

Fresh peripheral blood samples (5 mL) were obtained from 
antecubital vein in evacuated EDTA vials (X'acutainer, BD 
Biosciences, USA). DNA extraction was carried out using DNA 
extraction kit (Real Genomics, Banqiao City, Taiwan) from the 
whole blood. The quality and quantity of DNA were checked 
spectrophotometericaUy (Specord 200 spectrophotometer, Analy- 
tikjena, Germany). 



1.4 Confirmation of HIV Negative Status in ESNs 

Absences of pro-viral DNA of HIV-l was confirmed in each 
ESN individual using an in-house developed nested-PCR, as 
described previously [16]. 

1.5 Confirmation of Absence of CCR5A32 Mutation 

The most widely known genetic association with 32 bp deletion 
in the CCR5 gene [CCR5A32), which confers nearly complete 
protection from HIV-l infection in homozygous individuals was 
ruled out in ESNs by PGR as described previously [17]. 

1.6 Polymorphism at -238 (rs 361525) and -308 (rs 
1800629) Positions in TNF-a. Gene 

This polymorphism consists of a G to A substitution at positions 
-238 and -308 in the proximal promoter of the TNF-a. gene. For 
genotyping nested PGR was carried out resulting in amplification 
of 266 bp product (spanning -372 bp to -106 bp). The reaction 
mixture (20 ^lL) comprised of 50-100 ng of template DNA, 10 pM 
of each primer (Sigma Genosys, USA), 200 |J,M dNTP mix 
(Fermentas, Vilnius, Lithuania) and 1 unit of Taq DNA 
polymerase (NEB, USA) along with 2 \xL of lOX buffer (consisting 
of 0.1 M Tris-HCl, pH 8.8, 1.5 mM MgClg, 0.5 M KCl and 1% 
Triton X-100). The primers [18] TNF-Pl F: 5' -GAA GGA AAC 
AGA CCA GAG AC-3' and TNF-Pl R: 5' -ATC TGG AGG 
AAG GGG TAG TG-3' amphfied a 266 bp fragment. PGR 
conditions included initial denaturation at 94°G for 5 min 
followed by 35 cycles of denaturation at 94°C for 35 Sec, 
annealing at 57°C for 30 Sec, extension at 72°C for 30 Sec, with 
final extension for 5 min at 72°C. PGR was carried on a Master 
cycler gradient machine (Eppendorf, Hamburg, Germany). The 
PGR products were monitored by gel electrophoresis in a 2% 
agarose gel (containing 0.5 ng/mL ethidium bromide). DNA 
bands were visualized under UV light and digitally photographed 
using a gel documentation system (Syngene, Cambridge, England). 
Amplified products of first (external) PGR were used as template in 
the second (internal or nested) PGR to amplify 118 bp products. 
Nested PGR was carried out by using primers: TNF-Nc F: 5'- 
AGG GAA TAG GTT TTC AGG TCC ATG-3' and TNF-Bg R: 
5'-CAC ACT CCC CAT CCT CCC AGA TC-3'. The PGR 
conditions were same as mentioned above. Four micrograms of 
nested PGR product (118 bp) was digested at 37°C for 8 hours in a 
10 |J,L reaction volume containing 2.5 units of restri[:ti()n enzyme 
Ncol (NEB, USA) to detect -308G/A polymorphism [18] and 2.5 
units of restriction enzyme BgKl (NEB, USA) to detect -238G/A 
polymorphism, and then analyzed on a 15% non-denaturing 
polyacrylamide gel (PAGE). Restriction fragment length polymor- 
phism (RFLP) analysis for -308G/A SNP resulted in two bands (94 
and 24 bp) for GG genotype, three bands (1 18, 94 and 24 bp) for 
GA heterozygous and single uncut band (118 bp) corresponding to 
AA genotype. Similarly RFLP analysis for -238G/A SNP resulted 
in two bands (95 and 23 bp) for AA genotype, three bands (1 18, 95 
and 23 bp) for GA heterozygous and single uncut band (118 bp) 
corresponding to GG genotype. 

1.7 Polymorphism at -863 Position in TNF-a. Gene (rs 
1 800630) 

This polymorphism consists of a C to A substitution at position - 
863 in the proximal promoter of the TNF-a gene. The subjects 
were genotyped for this polymorphism by PGR-RFLP as 
described previously [19]. The region surrounding the polymor- 
phism was amplified using forward primer F:5'-GGG TCT GAG 
GAA TGG GTT AC-3' and reverse primer R:5'-CTA CAT 
GGG CCT GTC TTC GTT ACG-3'. Amplification was 
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performed in a reaction volume of 20 JtL containing 1 mM 
MgGl2, 0.2 mM of each nucleotide (Fermentas, Vilnius, Lithua- 
nia), and 1 U of Taq DNA polymerase (NEB, USA). PGR 
conditions included initial denaturation at 94°C for 5 min 
followed by 35 cycles of denaturation at 94°C for 30 Sec, 
annealing at 62°C for 1 min and extension at 72°C for 30 Sec, 
the final extension being at 72°C for 5 min. Four micrograms of 
amplified product (126 bp) was then digested with BsaAI 
restriction enzyme (NEB, USA) at 37°C for 8 hours, electropho- 
resed on 15% PAGE at 60 V and stained with ethidium bromide. 
The 1 26 bp band corresponded to the G wild-type C allele and a 
set of 103 bp and 23 bp bands corresponded to the variant A 
allele. 

1.8 Plasma TNF-a Level Estimation 

Plasma TNF-a concentration was estimated by Human GBA 
Thl/Th2 cytokine kit (Gytometric Bead Array, BD Biosciences, 
San Jose, GA, USA) according to the manufacturer's instructions. 
The lower limit of TNF-a detection using this kit was 2.8 pg/mL. 
Standard curves generated and plasma TNF-a concentrations 
were calculated using the BD CBA Software (BD Bioscience, San 
Jose, CA, USA). 

1 .9 Assessment of Blood CD4+ T-cell Count and Plasma 
Viral Load 

Absolute CD4 T-lymphocyte counts were estimated in the 

whole blood by flow cytometry using BD Tritest GD3 FITG/GD4 
PE/CD45 PerGP with BD Trucount tubes (BD Biosciences, San 
Jose, GA, USA) and plasma viral load (pVL) was estimated using 
GOBAS Amplicor HIV-1 monitor test, version 1.5 (Roche, NJ, 
USA) as per the manufacturer's guidelines. 

1.10 Mitochondrial Membrane Potential (A\|/m) 
Estimation 

Mitochondrial membrane potential (Av|/m) was measured flow 
cytometrically using 5, 5', 6, 6'-tetraethylbenzimidazolocarbocya- 
nine iodide JC-l) dye (Sigma-Aldrich, USA) [20]. Briefly, 2 mL of 
EDTA anti-coaggulated whole blood was treated with RBG lysis 
bulfer (Himedia, Mumbai, India) and leukocytes were washed with 
phosphate buffer saline (PBS, Himedia, Mumbai, India). GeU 
count was adjusted to 1 xlO^ ceUs/mL before staining with JG-1 
dye (final concentration 2.5 (xM) by incubation at 37°G (5% CO2 
humidified atmosphere) for 1 5 minutes in dark. Stained cells were 
washed with PBS and acquired immediately on flow cytometer 
(FACSGanto II). The data was analysed using software FAGSDiva 
Version 6.1.3 (BD Biosciences, USA). As a positive control for 
maximum depolarization of mitochondrial membranes, the 
mitochondrial uncoupling agent protonophore carbonyl cyanide 
p-(trifluoromethoxy) phenylhydrazone (FGGP) (Sigma-Aldrich, 
USA) was used at final concentration of 2 \lM. A total of 50,000 
events were acquired for each assayed sample. Lymphocytes were 
gated on forward and side scatter to exclude debris and non- 
lymphoid cells. JG-1 fluorescence was analysed on FLl and FL2 
channels for the detection of dye monomers (Green) and J- 
aggregates (Red), respectively. The ratio of red/green fluorescence 
reflected mitochondrial trans-membrane potential (A\|/m). The 
treatment with the protonophore FGGP resulted in decreased JC-1 
fluorescence ratio and served as a positive control for disruption of 
mitochondrial membrane causing change' in membrane' potential. 
The Av|/m is expressed as percent median fluorescence intensity 
(FL2:FL1/FL2 FCCP: FLl FGCP)xl00. 



1.11 Statistical Analysis 

Statistical analysis was performed using GraphPad Prism 
version 5.0 (GraphPad Software, Inc.). Power analysis was 
performed using Qrianto software (version 1.0; http://hydra.usc. 
edu/gxe). Hardy- Weinberg equilibrium (HWE) was examined for 
each SNP by online HWE calculator (http://www.oege.org/ 
software/hwe-mr-calc.shtml). The frequencies of genotypes and 
alleles were compared by Fisher's exact test. Fisher's exact test was 
also used to analyze association between HIV disease progression 
and TNF-ac variants. Odds ratio (OR) and its 99% confidence 
interval (G.I.) were calculated to assess the risk conferred by a 
particular allele and genotype. Haplotype analysis was performed 
using PHASE (http://www.stat.washington.edu/stephens/software. 
html) v2.1. Slope of CD4 T cell decline was calculated by linear 
regression model. Discrete and continuous variables were com- 
pared between cases and control subjects using unpaired Student's 
t-test. Two way ANOVA and Mann-Whitney test as appropriate. 
The Bonferroni correction for multiple testing was applied when 
necessary. A two sided p-value of <0.05 was considered 
statistically significant. 

Results 

A total of 260 healthy volunteers (HIV-1 seronegative, HGs), 
200 HIV-1 seropositive patients and 50 ESNs were genotyped for 
-863G/A, -308G/A, and 238G/A polymorphisms in the promoter 
region of TNF-a gene. Genotype distribution and allelic frequency 
were compared between HGs, SPs, FPs and ESNs. The 
distribution of the genotype frequencies at these loci were found 
to be in Hardy-Weinberg equilibrium. 

3.1 Characteristics of Study Population 

All individuals (5 1 0) included in this study were ethnically from 
the North Indian states. The demographic and clinical profile of 
individuals in dilferent study groups is summarized in Table 1 . No 
significant difference was found in the age and sex distribution 
between these study groups. AU ESNs were negative for HIV-l 
infection by serology as well as by DNA PGR and none of them 
was positive for the presence of CCR5A32 mutation (data not 
shown). For the HIV-1 seropositive individuals included in this 
study, we had the baseline values of viral load and GD4 count both 
at the time of registration in the ART clinic as well as at the time of 
recruitment in this study. At the time of registration the FPs had 
initially higher median CD4 count than SPs (617 vs 583 cells/nL, 
Range given in Table 1), but within one year there was a steep 
decline in GD4 count in FP group (slope of decline for GD4 count 
>-200 cells/|xL/Year, mean -464 cells/lxL/Year, Range -201 to 
-506 cells/ H-L/Year), whereas in SP group the rate of decline was 
only marginal (slope of decline <-32 celIs/(i.L/Year, mean -18 
ceUs/nL/Year, range 20 to -32 cells/ jtL/Year) even after more 
than three years of HIV detection. Significant difference was 
found in the slope of GD4 count decline in FPs as compared to SPs 
(Unpaired T-test, jft<0.000 1). The HIV plasma mean log viral load 
was also significantly (Unpaired T-test, /)<0.0001) higher in FPs 
(5.4 Logio RNA copies/mL) as compared to SPs (3.4 Logio RNA 
copies/mL, range given in Table 1). 

3.2 TNF-a. Gene -238G/A Polymorphism 

On genotyping of -238G/A polymorphism (Figure lA) in the 
promoter region of TNF-a. gene we found genotype distribution of 
GG (91%) and GA (9%) in the healthy controls with allele 
frequency of G allele at 95.38% and A allele at 4.62%, which was 
considered as wild type distribution (Table 2). As compared to this 
no significant difference was found in genotype and allele 
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Table 1. Demographic and clinical characteristics of the population analyzed. 





Characteristic 


HCs (n = 260) 


SPs (n = 100) 


FPs (n = 100) 


ESNs (n = 50) 


Sex (M/F ratio) 


1.05 


1.2^ 


1.1" 


0.87' 


Age (years) {median, range) 


27 (19-48) 


34.5 (18-52) 


35 (21-51) 


33 (24-42) 


Duration of HIV-1 infection (months) (mean) 




60 (38-96) 


12 (8-14) 




CD4 T-cell count (cells/^il) (median and range) at the time of HIV Diagnosis 




583 (419-1116) 


617 (417-733) 




CD4 T-cell count (cells/^il) (median and range) at the time of recruitment 


862.5 (463-1471) 


495 (411-884)11 


202 (142-227)11 


600 (465-884) 


Plasma HIV-1 RNA (copies/ml) (mean Log Value and range) 




3.4 (2.7-4.1)* 


5.4 (4.87-5.85)* 




Exposure to HIV-1 History 


Sexual 




82 (%) 


88 (%) 




Blood Transfusion 




4 (%) 


2 (%) 




Intravenous Drug User 




2 (%) 


3 (%) 




Other 




1 2 (%) 


7 (%) 





a, SPs vs HCs, p>0.C5; b, FPs vs HCs, p>0.05; c, ESNs vs HCs, p>0.05. 
*Viral Load, FPs vs SPs, p<0.00C1. 
IICD4 Count, FPs vs SPs, p<0.0001. 
doi:1 0.1 371 /journal.pone.0098020.t001 



frequency of this polymorphism in HIV- 1 seropositive individuals. 
Similarly there was variation in distribution but difference was not 
significant within HIV seropositive groups viz. SPs vs FPs or ESNs 
group (Table 3). 

3.3 TNF-a Gene -308G/A Polymorphism 

The genotype distribution of GG (90%), GA (9%) and AA (1%) 
was observed in healthy controls with allele frequency for G allele 
as 92.42% and for A allele 5.58%. No significant difference was 
found in genotype and allele frequency when compared with 
overall HIV-1 seropositive group (Table 2). However, the 
frequency of GA genotype, which is known to be associated with 
higher production of TNF-a cytokine, was significantiy higher in 
fast progressors (FPs) as compared to slow progressor (SPs) group 
(/)<0.01, OR= 3.43). Further comparison between allele frequen- 
cies between FPs vs SPs also showed A allele to be significantly 
over-represented in FPs (p<O.Ol, OR=3.39) (Table 3). RFLP 
analysis observed on genotyping TJVF-oc -308G/A is shown in 
Figure IB. 

3.4 TNF-a Gene -863 C/A Polymorphism 

Genotype analysis of TMF-cc -863C/A SNP revealed, CC (50%), 
CA (43.85%) and AA (6.15%) in healthy controls with frequency 
of C allele as 71.92% and of A allele as 28.08%. No significant 
difference was observed in the frequency distribution of these 
alleles among HFV positive groups when compared to HC 
(Table 3, Figure IC). 

3.5 Haplotype Association with Disease Progression 

In order to evaluate the combined effects of the three SNPs in 
the TMF'CC gene, haplotype analysis was conducted using PHASE 
software. Of the 8 haplotypes observed in the patients and controls 
recruited for the study, only three haplotypes: CGG, CAG and 
AGG were represented at frequencies >5%. These accounted for 
more than 99% of the total haplotypes. Among these, the 
frequency of wild type haplotype CGG in healthy controls (65%)) 
was not significantly different from that in SPs (68%). However, 
the CAG haplotype, which is high producer of TNF-a, was 
significantly overrepresented in FPs {p<0.03, OR= 1.93, 99% CI 
(0.88-4.21)] as compared to HCs and SPs \p<Om, OR = 3.00, 
99% CI (0.98-9.20)] (Table 4). The frequency of low producer 



haplotype AGG however, was not significantly different among 
various study groups. 

3.6 Plasma TNF-a Concentrations 

TNF-a levels were measured in the plasma samples of 18 
subjects in each group. The plasma TNF-a levels were signifi- 
cantiy (/)<0.001) higher in FPs (20. 9± 1.4 pg/mL) as compared to 
SPs(7.9±l.l pg/mL) and ESNs (4.1±0.6 pg/mL,^<0.0001). We 
further investigated the effect of haplotypes on TNF-a levels within 
FP group. Individuals having CAG (high producer) haplotype had 
significantiy (/)<0.01) higher levels (Mean, 27.36±2.4 pg/mL) of 
TNF-a as compared to carriers of low producer AGG (17 ±1.5, 
pg/mL) & wild type CGG haplotype (18.26±2.4 pg/mL, /)<0.01) 
(Figure 2). However, the comparison between haplotypes within 
SP group showed a significant difference between high producer 
and low producer haplotypes only (CAG vs AGG: 9.98±1.6 pg/ 
mL vs 6.08± 1.5 pg/mL, /;<0.01). Further, in FPs, positive 
correlation was also observed between plasma TNF-a concentra- 
tion and plasma viral load (Pearson correlation, r = 0.77, p< 
0.0002). 

3.7 TNF-a Genotypes and Clinical Parameters 

We tried to investigate the distribution of TNF-a. haplotypes and 
baseline viral load and CD4 T-cell counts within FP group. High 
TNF-a producing haplotype CAG within FP group had signifi- 
cantly higher mean log viral load as compared to carriers of wild 
type haplotype (CGG) (Mean, 5.8 vs 5.4 Logio RNA copies/mL, 
Unpaired T-test, /)<0.001) as well as low producer AGG and 
faster decline of CD4 counts as observed by steeper decline of 
CD4 slope (mean slope of decline for CD4 count -484 cells/ |a.L/ 
Year vs -464 cells/nL/Year) (Unpaired T-test, /)<0.01) (Figure 3 
A, B). 

3.8 Haplotypes and IVlltochondrial IVlembrane Potential 

Since, the high producer genotypes and haplotypes combined 
with high plasma TNF-a concentration were found to be 
associated with faster decline in CD4 count, we investigated 
further, whether these genotypic differences could influence the 
apoptosis of lymphocytes in terms of decrease in mitochondrial 
membrane potential (Av|/m) in HIV infected individuals (Figure 4). 
In order to address that, subjects (n = 5) were recruited in each 
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118 bp 
95 bp 




118 bp 
94 bp 

20 bp ladder 



126 bp 
103 bp 

50 bp ladder 



Figure 1. Restriction fragment lengtli polymorpiiism (RFLP) analysis of PCR products. (A) Representative ethidium bromide stained (15%) 
non-denaturing polyacrylamide gel (PAGE) sliowing amplified gene product of TNF-x (-238G/A) 118 bp uncut and RFLP pattern observed after 
digestion with 6g/ll. (Lane M -20 bp ladder. Lane 2 and 8-GA genotype. Lane 3 to 7-GG Genotype. (B) Representative non-denaturing PAGE (15%) 
showing amplified gene segment of TNF-a (-308G/A) 118 bp uncut and RFLP pattern observed after digestion with Nco\. Lane 1-uncut 118 bp 
product. Lane 2 and 6-GG Genotype, Lane 3 and 5-GA genotype, Lane 4-AA Genotype, Lane M-20 bp Ladder. (C) Representative non-denaturing 
PAGE (15%) showing amplified gene segment of 7WF-a (-863C/A) 126 bp uncut and RFLP pattern observed after digestion with fisoAI. Lane 1- 
Negative Control, Lane 2-uncut 126 bp. Lane 3-AA Genotype, Lane 4-CA Genotype, Lane 5, 8 and 9-CC Genotype, Lane 6 and 7-CA Genotype, Lane M- 
50 bp ladder. 

doi:1 0.1 371/journal.pone.0098020.g001 



group based on the haplotypes [GAG (high TNF-a producer) vs 
CGG (low TNF-a producer)]. Significandy low mean A\|/m was 
observed in HIV-1 infected patients as compared to healthy 
controls (669±34 vs 920±45, Unpaired T-test, /)<0.003). Also, 
Av|/m was significandy lower in FP group as compared to SP group 
(493±25 vs 703±24, Unpaired T-test /;<0.002). Further, when we 
compared A\|/m of different haplotypes within each group or 
between groups (Figure 5), we found that individuals having GAG 
haplot\pe had significantly lower mean A\|/m as compared to 
carriers of wild t^pe (CGG) haplotype within FP group (41 7±22 vs 
571 ±28, Two way ANOVA, p<O.Ol), indicating higher cell death 
among high TNF-a producing haplotype carriers as compared to 
wild type. Similarly in case of SP group also, individuals having 
GAG haplotype had lower A\|/m as compared to carriers of wild 
type (CGG) (638±19 vs 768±29, Two way ANOVA, /><0.05). 

Discussion 

The present study evaluated the frequency of TMF-a promoter 
polymorphisms and their association with HIV-1 disease progres- 
sion or risk of HIV infection, plasma levels of TNF-a and 
mitochondrial membrane potential in ethnically defined North 
Indian population. The TNF-a gene is located within the class III 
region of the major histocompatibility complex (MHC) between 



HLA-B and DR, and its expression in the host is controlled in a 
regulated manner at the transcriptional as well as post transcrip- 
tional level. Polymorphism in the TJVF-a promoter region has been 
shown to play an important role in the development of 
autoimmune diseases and several infections like HBV, HCV, 
leishmaniasis, etc. [5] [21] [22,23]. The -308A aUele in the TNF-a 
promoter is a G/A transition and is the most studied SNP in this 
gene. Our data indicate that TNF-a -308 G/A polymorphism is 
associated with faster progression of HIV disease both at genotype 
and allele frequency level. Although, when Bonferroni's correc- 
tions are applied for multiple comparisons, the significance level 
comes at border line only(corrected /; values: p<O.Ol, OR = 3.43, 
p<O.Ol, OR=3.39). Haplotype analyses also suggest that the 
combination of some genetic variants within the TNF-a gene (e.g. 
GAG haplotype), conforming to higher production of the cytokine 
was represented in significandy higher number of FPs as compared 
to SPs, suggesting that the high producers of TNF-a would show 
faster disease progression. These results do agree with some 
previous reports, which have shown an association between TNF-a 
SNPs -238G/A and -308G/A and the risk of HIV-1 infection [13] 
[14] while not with others which observed no significant 
association [15]. This discrepancy may be attributed to the 
different ethnic background of assessed populations in these 
studies. 
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Furthermore, we observed a significantly lower mean lympho- 
cyte mitochondrial membrane potential (Av|/m, an indicator of cell 
apoptosis) in FPs as compared to SPs and even within group (FPs) 
carriers of CAG haplotype as compared to carriers of wild type 
haplotype (CGG), indicating that the individuals having high 
TNF-a producing genetic background wall have faster disease 
progression due to predisposition to higher rate of cell death. 

In the current study, we have found almost similar TNF-a. 
-238G/A and -308G/A allele frequencies as reported earlier from 
North Indian population in studies related to Mycobacterium, 
tuberculosis infection and Type-1 Diabetes [24] [25]. It seems 
therefore that -238A is a rare allele in Indian population (with 
frequency varying between 2 to 5%) as compared to European 
population which has been reported to be more than 8% [19] . It is 
also important to note that the population frequency of TNF-a. 
-308A in the Caucasian population is around 17-22%, which is 
higher than that observed in the Asian populations (5-6%), which 
is rare as compared to -308G allele (found in 93%), suggesting 
large variation in the allele frequencies among different ethnic 
groups [26]. We did not find any significant difference between the 
frequencies of any of the TNF-a alleles or haplotypes in the ESNs 
and HIV infected groups, indicating no association of these 
polymorphisms with risk of HIV infection in North Indian cohorts. 
In contrast to our results, a recent study from Western India 
(Mumbai) showed a higher frequency of -238GG genotype in 
exposed uninfected (EU) individuals as compared to those who 
acquired infection, indicating a possible association with HIV 
transmission [27]. Another study on Spaniard cohort showed no 
association with any of the individual polymorphisms at -238(G/ 
A), -308(G/A) and -863(C/A) loci with the vulnerability to HIV 
infection, however the haplotype analysis in their study indicated 
that haplotype GAG was more frequently represented in the EU, 
suggesting the combination of polymorphisms within TNF-a gene 
may positively modulate the risk of HIV infection in an individual 
[19]. The discrepancy in findings may be due to the ethinic 
differences in the studied populations. So, further studies at 
genotype and haplotype levels are needed for a logical explanation 
regarding association of TNF-a polymorphism with risk of 
infection. 

Similar observations have been made in relation to different 
infections, like in case of HBV infection, the presence of -308A and 
-863G was found to be associated with viral clearance [5], whereas 
SNP -238G has been shown to be associated with protection 
against chronic hepatitis in the Chinese and German population. 
In hepatitis C virus infection, the haplotype -863C/ — 308G was 
associated with viral persistence in African Americans. Similarly 
the existence of -308A SNP in pediatric population positively 
correlated with the severity of dengue viral infection in children 
[28]. With regard to HIV-1 infection, studies have shown 
association of TNF-a -308A allele with HIV-1 -associated dementia 
and other HIV therapy related complications like Lipodystrophy 
Syndrome in infected patients [10]. 

The immuno-regulatory response of the host influences the 
pathogenesis of HIV-1 infection, triggering monocytes, macro- 
phages, and natural Idller cells to produce TNF-ot [29]. We have 
found significantly high TNF-Gt levels in plasma of FPs as 
compared to SPs and ESNs and there is a positive correlation 
between HIV-1 viremia and TNF-a levels in plasma of HIV-1 
infected individuals. This suggests that reducing TNF-a levels may 
also help in reducing the HIV- 1 viral load. In excess, TNF-a may 
cause severe inflammatory damage and toxicity, making control of 
its production and secretion highly important. Regulating its 
release serves as a potential means of therapy for HIV- 1 and other 
diseases. TNF-a can also induce other pro-inflammatory cytokines 
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Figure 2. Plasma concentrations of TNF-a in various study groups. Plasma concentrations of TNF-ot (Mean±SEIVI) in FPs, SPs and ESNs. High 
TNF-a producing haplotype CAG (dark bar), medium TNF-a producer haplotype CGG (light gray bar) and low TNF-a producing haplotype AGG (gray 
bar). 

doi:10.1371/journal.pone.0098020.g002 
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Figure 3. Plasma viral load and CD4 counts in different TNF-a haplotypes (CGG, CAG and AGG). (A) Plasma Log^o RNA copies/mL 
(Mean±Range) in TNF-a haplotypes, CGG (dark bar), CAG (light gray bar) and AGG (gray bar). (B) CD4 counts at different time intervals in TNF-y. 
haplotypes, CGG (dark bar), CAG (light gray bar) and AGG (gray bar). Values shown as (Median±SD). Significance levels: **p<0.01; *p<0.05. 
doi:1 0.1 371/journal.pone.0098020.g003 
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Figure 4. Gating strategy for the analysis of JC-1 stained PBIVICs. Lymphocytes were gated according to morphological parameters (not 
shown in figure), JC-1 aggregates and monomers were analyzed on FL2 (PE) and FL1 (FITC) channel respectively, (a) showing JC-1 stained 
lymphocytes of HCs (P2 gate, 96.4%), (b) protonophore FCCP treated lymphocytes as a positive control showing 84.2% cells having reduced 
mitochondrial membrane potential, (c) and (d) JC-1 stained lymphocytes in a representative samples from individuals in 'fast progressors' group 
having CGG and CAG haplotypes showing 26% and 37% lymphocytes dying respectively. Change in mitochondrial membrane potential (Av|/m) was 
expressed as percentage of median fluorescence intensity (FL2: FL1/FL2 FCCP: FL1 FCCP) xlOO. 
doi:1 0.1 371 /journal.pone.0098020.g004 
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Figure 5. Median fluorescence intensities of total lymphocyte mitochondrial membrane potential (Ai|/m). Bar diagram showing median 
fluorescence intensity (MFI) of total lymphocyte mitochondrial membrane potential (Av|)m) in FPs and SPs having CAG (n = 5) (dark bar) and CGG 
(n = 5) (light gray bar) haplotypes (shown as IVlean±SEM). 
doi:1 0.1 371 /journal.pone.0098020.g005 



such as IL-6 and IL-8, which aid in the up-regulation of viral 
rephcation [30] [31]. TNF-a is secreted during the early phase of 
acute inflammatory diseases. Its pathogenic role in HIV-1 
infection involves activation of nuclear factor kB (NF-kB), 
stimulating viral replication as well as apoptosis of T lymphocytes, 
both the factors contributing to faster disease progression. Tissue 
and plasma samples of hosts express high levels of TNF-ot, 
contributing to fever, anorexia, and other symptoms of HIV/ 
AIDS. These observations suggest that TNF-a must be targeted at 
an appropriate time during the disease to prevent its progression to 
the chronic stage. Local effect of the cytokine may be beneficial to 
the host, so monitoring its development is critical. 

The current study is novel in terms of our study population as 
the three highly unique HIV- 1 subtype-C infected cohorts from a 
previously not seriously investigated ethnic population, were 
studied for the expression of polymorphic alleles that could serve 
as clinical biomarkers for the progression of HIV-1 infection. To 
the best of our knowledge this is the first study from North Indian 
population on TJVF-oc promoter gene polymorphism in HIV-1 
infection having different clinical forms. Also, we have tried to find 
association at the haplotype level as haplotype-based approaches 
are more powerful than single SNP approaches in the setting of 
multiple susceptibility loci. However, additional multi-centric 
studies with larger number of patients from similar cohorts are 
warranted to establish the detailed mechanisms underlying the 



functional effects of these genetic polymorphisms. Also, HLA 
typing of these HIV- 1 infected patients is required as recently in 
genome wide association studies (GWAS) it has been shown that 
the TNF-a. SNPs are in linkage disequilibrium with MHC class I 
gene. It becomes more important in the light of well documented 
association of some HLA types with the rate of HIV disease 
progression [32,33]. None-the-less the data given in this study 
provides significant information regarding genetic association with 
progression of disease, which might be very helpful in designing 
control strategies in such populations. 
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